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Precisely measured electron-proton elastic scattering cross sections [Pliys. Rev. Lett. 105,242002 
(2010)] are reanalyzed to evaluate their strength for determining the rms charge radius {Re) of the 
proton. More than half of the cross sections at lowest are fit using two single-parameter form- 
factor models, with the first based on a dipole parametrization, and the second on a linear fit to 
a conformal-mapping variable. These low-Q^ fits extrapolate the slope of the form factor to Q^=0 
and determine Re values of approximately 0.84 and 0.89 fm, respectively. Fits spanning all Q^, 
in which the single constants are replaced with cubic splines at larger Q^, lead to similar results 
for Re- We conclude that the scattering data is consistent with Re ranging from at least 0.84 to 
0.89 fm, and therefore cannot resolve the discrepancy between determinations of Re made using 
muonic and electronic hydrogen-atom spectroscopy. 

PACS numbers: \pacs{06.20. Jr, 13.40. Gp, 14.20 . Dh, 25.30. Bf > 


Recent measurements of the n=2 energy intervals of 
muonic hydrogen, when compared to precise QED theory 
for this exotic atom, lead to a determination mm of the 
rms charge radius of the proton (Re) of 0.84087(39) fm. 
This value disagrees by 4.5 standard deviations with a 
value of 0.8758(77) fm obtained from a similar compari¬ 
son [2 between QED theory and several precision mea¬ 
surements in the ordinary hydrogen atom. A third de¬ 
termination of i?E can be obtained from precise measure¬ 
ments of the the cross sections for elastic scattering be¬ 
tween electrons and protons. The most precise electron- 
proton scattering experiment is the recent measurement 
[1] of the MAMI collaboration, and their analysis [S] leads 
to Re= 0.879(8) fm, which disagrees with the muonic hy¬ 
drogen value by 4.6 standard deviations. 

CODATA [3] uses a combination of the scattering and 
hydrogen values to obtain Re, and its value differs from 
the muonic hydrogen value by 7 standard deviations. 
This disagreement has now widely been referred to as 
the proton size puzzle [S]. Many papers have discussed 
this puzzle, including many that have proposed physics 
beyond the standard model [?]• 

Because of the importance of the electron-proton scat¬ 
tering data to this puzzle, the data has been extensively 
scrutinized and discussed 13 [HHin]. The present work 
reanalyzes the scattering data and concludes that it is 
consistent with a much larger range of Re values than 
obtained by others. This range makes it consistent with 
both the hydrogen and muonic hydrogen determinations 
of Re, therefore removing one component of the proton 
radius puzzle. 

Our analysis uses two separate simple one-parameter 
form-factor models. Both of these simple models fit well 
to the low-Q^ data, but the two give discrepant valnes 
for i?E- Since neither model can be ruled out, the un¬ 
certainty in Re must, at minimum, be expanded to en- 
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compass both values. Generalizations of both models fit 
well to the entire MAMI data set, and give similarly dis¬ 
crepant values for Re- 

The differential cross section for elastic scattering of 
an electron of energy E scattering by an angle 6 from 
a stationary proton, (after taking into account radiative 
corrections and two-photon exchange) can be written |21] 
in terms of the squares of the electric and magnetic form 
factors (Ge(Q^) and Gm(Q^)): 


(Tred = (1 + r) 


da 

dD 


dD 



( 1 ) 


where dcTMott/dfl is the Mott differential cross sec¬ 
tion, e=(l -I- 2(1 -I- T)tan^|)“^, T=—t/{dm'^), and 
t=—Q^={pi —Pf)^, with Pi and pf being the initial and fi¬ 
nal four-momenta of the electron. Here, nip is the proton 
mass, and we use units with h = c = 1. 

In principle, the quantity of interest for this work. 


Re = 




1 ry red 

dt 


V dt 


t=o 4m^ ’ 


( 2 ) 


could be determined using Eq. Q from sufhciently- 
precise measurements of dcr/dH for small values of Q^. In 
practice, for the existing set of measurements, an extrap¬ 
olation to Q^=0 is required, and, for this extrapolation, 
a functional form for G| and G^ of Eq. ([^ must be 
assumed. 

The dipole form of the form factor has been used for 
many decades [22] and it approximates Ge and Gm/Pp 
( where Pp is the magnetic moment of the proton in units 
of nuclear magnetons) as: 



A second approximation for the form factors is to use a 
Taylor expansion in t about t=0. This Taylor expansion 
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has a limited radius of convergence due to a negative- 
pole at t=Am^ (where is the mass of the charged 
pion), which results from the two-pion production thresh¬ 
old. A conformal mapping variable |23| : 


y/tc ~ t — y/i^ 
y/tc"— t y/i^ 


(4) 


with tc=4m^ leads to a much larger radius of conver¬ 
gence. Thus, 


G| = 1 - cez, and Gm/a^p = 1 - cm^ (5) 


are good approximations to the form factors at low Q^. 

Other functional forms for Ge and Gm have been used, 
to extrapolate to <5^=0 to determine Re- These other 
forms include polynomials in t ms], polynomials in z 
[13 na 1181. inverse polynomials in t 00, dipole func¬ 
tions (Eq. Q) times polynomials in t [5], dipole functions 
plus polynomials in t [3, cubic splines in t 00, dipole 
functions times cubic splines in t [3, continued fractions 
in t |10j , and the Friedrich-Walcher parametrization (two 
dipole functions plus two symmetric gaussian features) 
[ 3 . In this work, we restrict ourselves to the forms of 
Eqs. ^ and ([^ for low-Q^ data, and extensions of these 
forms for higher-Q^ data. 

The highest-accuracy e-p scattering experiment 0by 
the MAMI collaboration yields 1422 cross sections (with 
typical relative uncertainties of 0.35%) spanning a range 
of 180 MeV < E <855 MeV and 16° < 6» < 135.5°, 
corresponding to 0.0038 GeV^ < < 1 GeV^ and 

0.06 < e < 1. The 1422 cross sections are divided into 
34 data groups, with each data group having a separate 
normalization constant. These normalization constants 
are known to an absolute accuracy of a few percent, and 
are related to one another in such a way that there are 
only 31 independent constants 00. 

The normalization constants add a further complica¬ 
tion to the Q^=0 extrapolation needed to determine Re- 
The few-percent absolute accuracy of the measured cross 
sections is not sufficient for performing a precise extrap¬ 
olation, and thus the 31 normalization constants need to 
be floated when performing least-squares fits of the entire 
data set for this extrapolation. We include these normal¬ 
ization constants in all of our fits. Floating these con¬ 
stants adds considerable flexibility to the extrapolations. 
Although we do not impose the few-percent absolute un¬ 
certainty of the normalization constants in our fits, all of 
our fits return constants near unity and well within this 
few-percent uncertainty. 

Other least-squares fits H13 m USUI [13 of this data 
use seven- to twelve-parameter models for Ge and for 
Gm, and obtain least-squares fits with reduced values 
of as low as 1.14 for fitting the 1422 data points. The 
1.14 value is much too large for the number of degrees of 
freedom in the fit, but can easily be explained by either a 
7% underestimation of the uncertainties, or a systematic 
effect that is not fully accounted for. We only include fits 
that have a reduced x^<1.14 in this work. 


As indicated by Eq. ([^, the rms charge radius of the 
proton is a small-concept. Thus, if possible, it should 
be determined from low-Q^ data. Therefore, we attempt 
to make a determination of Re using fits to only the 
lower-Q^ data. In addition to the fact that such fits use 
data nearer to Q^=0, the fits have the advantages that 
simpler, fewer-parameter models can be used for G|. and 
G^, and that, since fewer of the data groups are used, 
fewer normalization constants need to be included in the 
fits. 

Fits using single-parameter models for the form fac¬ 
tors are shown in Fig. These fits include data with 
Q^<Qmax=0-l GeV^, from 19 data groups, which re¬ 
quire 17 normalization constants. A data group is only 
included if there are more than 10 data points in the 
group, and a total of 761 of the 1422 cross sections (53%) 
are used. The value of Re is directly obtainable from the 
slope of the curves in Fig. [l]at Q2 _o, and the fits provide 
the necessary extrapolation to Q^=0- 

Fig. [^a) shows a fit using Eq. Q and the one- 
parameter dipole form factors of Eq. (§. The 
reduced for the fit is 1.11, and the fit re¬ 
turns i?E=(12/&_E)^/^=0.842(2) fm and i?M = (12/6M)^^^ 
=0.800(2) fm. 

A second fit to the same data uses the one- 
parameter linear model (in z) of Eq. ([^. This fit 
is shown in Fig. Sb). , It also has a reduced 
of 1.11 and gives i?E=(|cE/tc)^^^=0.888(l) fm and 
i?M=(|cMAc)^/"=0.874(2) fm. 

Figure shows (red bands labeled Eq. ([^ and Eq. ([^ 
at the left of the figure) the error bands for Re for the 
dipole and linear fits versus the cutoff Q^ax- The figure 
includes the range of Q^^x for which a reduced x^<1.14 
is obtained. 

The electric form factors predicted from the two fits of 
Fig. [3 are shown in Fig. [3 Also plotted in the figure are 
other low-Q^ measurements of these form factors (often 
referred to as the world data, as summarized in Ref. [24] 1. 
It is clear from this figure (and from the calculated 
for the comparison betwen the data and the two curves) 
that the form factor from either fit is also consistent with 
these other measurements. 

One concern that could be raised about the single¬ 
parameter fits, which are based only on data with < 
is that they may lead to inconsistencies for data 
with (5^>(5aiax' Since the low-Q^ fits of Fig. deter¬ 
mine 17 of the 31 normalization constants, and since 
data groups using these normalization constants include 
measured cross sections with > Q^ax, iFe fits have 
a direct impact on data not included when fitting. It 
could, therefore, be possible that data at the same value 
of > Qmax from two of these data groups could be 
made inconsistent when these normalization constants 
are used. 

To ensure that such inconsistencies do not take place, 
we extend the fits of the previous section to include all 
of the MAMI data. Such an extension also allows for a 
direct comparison of the quality of our fits to the quality 
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FIG. 1. (Color online) Fits to the low-Q^ data of Ref. [1] using single-parameter models for the form factors in Eq. In 
(a), the dipole model of Eq. ^ is used, and in (b) the linear model of Eq. is used. Separate plots are required since the 
fits return different normalization constants, and therefore the cross sections take on slightly (~1%) different values for the two 
plots. The i?=180 MeV data and fits are continued in the insets, and the other energies are offset for clarity of presentation. 


of the fits performed by others [H [51 HUl HU [151 US] who 
also include all of the MAMI data. The extended fits 
include a range of in which the functional form of 
and becomes more complicated, and, as with the fits 
performed by others, more parameters are necessary to 
obtain a good fit. 

The fit using Eq. ([^ can be generalized by allowing 
the constants 6 e and &m to become functions of z. We do 
this by using cubic splines, with the b values each being 
a constant for z < Zo=0-1, and with iVkn equally spaced 
knots between zq and Zmax- The values of &, and their 


QLx [GeV^] 

0.1 0.2 0.3 0.4 0.5 0.6 0.8 1 



FIG. 2. (Color online) The range of Re predicted from 
the dipole fits (lower bands) and conformal-mapping fits (up¬ 
per bands) as a function of the cutoff value Qmax- The red 
bands correspond to single-parameter dipole and linear fits 
of Eqs. and ([^, and the other colors show extensions to 
these fits for which the single parameters are replaced with 
cubic splines (with the number of nodes, Akn, as indicated) 
at larger . All fits shown have a reduced x^<1.14. 


first and second derivatives are continuous at Zq and at 
the other knots. The number of knots needed to achieve 
a good fit increases with increasing Zmax- 

The results of such fits with A^kn=l to 10 (2 to 11 pa¬ 
rameters per form factor) are given in Fig.[^ Again, only 
fits with a reduced x^<1.14 are shown. The fits return 
values of i?E of approximately 0.84 fm for all values of 
Qmaxj similar to the single-parameter dipole fit. The fits 
at the right of the plot include all of the MAMI data, 
and still return a value near 0.84 fm. 

Equation ([^ would predict negative values for G| and 
G^ at larger z. This problem can be avoided by using a 


Cf [GeV^l 



FIG. 3. (Color online) Electric form factor Ge determined 
from the single-parameter fits of Fig The blue curve is 
from the dipole-model fit, and the red curve is from the linear- 
model fit. Also shown on the plots are the low-Q^ Gb data 
from other experiments (from [24]), and these data are con¬ 
sistent with the Ge obtained from either fit. 
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denominator to cause a cutoff at higher z values: 


G 


2 _ 
E — 


1 - Ce^: 

1 - (cez)^ 


Gm 1 - cmz 


( 6 ) 


We use P=4 for our fits, but any P from 4 to 14 gives 
similar results. This function is very nearly linear up to 
2;=0.2, while avoiding negative values at larger z. A fit of 
the data to these form factors is performed by allowing 
ce and cm to become functions of z by using the same 
form of cubic splines as are used for 6e and &m- 

The results of such fits with A^kn=l to 10 are also given 
in Fig. The fits return values of Re of approximately 
0.89 fm, similar to the single-parameter linear fit. 

Fig. [^clearly indicates that the two types of fits pro¬ 
duce values of i?E that disagree with each other. Since 
either type of fit gives an extrapolation to zero that is 
equally valid, and since the quality of the fits are similar, 
either value of Re is possible. Therefore, at best, the 
determined value of Re can range from 0.84 to 0.89 fm. 
At worst, other valid extrapolations could lead to even a 
wider range of possible values, leading to an even larger 
range for Re- 

It is not the aim of this work to determine the rms 
magnetic radius of the proton, Rm, but we note that the 
values from our fits range from about 0.80 to 0.90 fm, and 
therefore this work cannot determine Rm to any better 
than this range. The consistency of the fits presented 
in this work can be checked by comparing the quantity 
Re + Rm to the prediction from hydrogenic spectroscopy. 
The hydrogen hyperfine interval determines to 

be 1.35(12) fm^, and the muonic hydrogen hyperfine in¬ 
terval determines i?| -I- to be 1.49(18) fm^ The 
two determinations are consistent, and their weighted av¬ 
erage gives i?| + = 1.39(10) fm^. The dipole fit of 

Fig.j^^gives R'^ + R'^ = 1.349(4), whereas the linear fit 
of Fig. mb) gives 1.553(4). The dipole fit is in excellent 
agreement with the spectroscopy result, and the linear 
fit shows only a mild 1.6 standard deviation discrepancy. 
Similar comparisons using the extended cubic-spline fits 
lead to a similar level of agreement. 

The extent to which two-photon exchange (TPE) af¬ 
fects the extraction of i?E from the MAMI data has been 
debated in the literature [aiilSlIIIlIIHlIMa- The cross 
sections given in Ref. [S] were corrected by the Coulomb 
corrections (Feshbach corrections [33]) in place of the full 
TPE corrections. In this work, these Coulomb correc¬ 


tions are removed and and replaced with TPE corrections 
calculated following the prescription of [341 [35] . This re¬ 
placement leads to correction factors of between 0.997 
and 1.003 for the data of Eig. and of between 0.978 
and 1.003 for the full MAMI set. The correction factors 
agree with those shown in Fig. 5 of Ref. [16] to within the 
0.03% accuracy readable from their figure. To test how 
sensitive our analysis is to TPE corrections, we repeat our 
full analysis using the low-Q^ TPE approximation of Ref. 
[36] and the Feshbach correction in place of the full TPE 
correction. Fig. [^ shows that using the Feshbach cor¬ 
rection would underestimate Re by only 0.004 fm, while 



FIG. 4. (Color online) The sensitivity of the fit results 
of Fig. [^ to two-photon exchange corrections is illustrated 
by redoing the analysis with poorer approximations for these 
corrections. The change in Re when the Feshbach correc¬ 
tions (solid curves) and a low-Q^ approximation [3^ (dashed 
curves) are applied shows that the radius extracted from both 
the dipole model (a) and conformal-mapping model (b) is not 
very sensitive to these corrections. 


the low-Q^ approximation would change Re by less than 
0.001 fm. We conclude that, although the best available 
TPE corrections should be used, the sensitivity to using 
poorer approximations is small in our analysis. 

In summary, we have reanalyzed electron-proton elas¬ 
tic scattering data using simple fits to the lowest-Q^ half 
of the data, and cubic-spline extensions of these fits at 
higher Q^. We find that the required extrapolation to 
Q^=0 can lead to values for the rms charge radius i?E 
ranging from 0.84 to 0.89 fm. This range does not re¬ 
solve the discrepancy between determinations of i?E from 
muonic hydrogen [Hi] and ordinary hydrogen [3]. 
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